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Abstract 

Schwinger operator method is applied for studying CP-odd pure gluonic effective 
Lagrangian in the Standard Model at three-loop level. The induced 0-term vanishes 
by the same reasons as EDMs of quark and W-boson to two-loop approximation. A 
simple way is found to demonstrate these cancellations. All other terms of the effective 
Lagrangian acquire non-vanishing contributions. The effective operator of dimension 
six, Weinberg operator, is calculated explicitly. The corresponding contribution to 
the EDM of neutron is much smaller than that comes from large distances. 
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1 Introduction 



The Kobayashi-Maskawa (KM) model looks now as the most natural description of CP- 
violation. It describes properly CP-odd phenomena in the decays of neutral if-mesons and 
predicts extremely tiny CP-odd effects in the flavour- conserving processes. Though its pre- 
dictions for the electric dipole moments (EDM) of elementary particles are far beyond the 
present experimental facilities, the corresponding theoretical investigations are of certain 
methodological interest. 

The subject of this work is the calculation of the CP-odd effective gluonic Lagrangian 
which appears in the Standard Model as a result of integration over quark and W-boson 
modes. This Lagrangian can be naturally expanded in the series of gluon field operators of 
increasing dimension: 

S eff = J d*x£ eff (x) = J d'x (c ig 2 G; u Gl + c 2 g 3 r bc G a a ,G b ^ a + ...) , (1) 

where g is the chromoelectric charge, G af3 = l/2G^ u e^ ua /3. It is assumed that the charac- 
teristic loop momenta are much larger than the inverse scale on which field fluctuations 
occur. The first term in (|1|) represents the induced (9-term, perturbative contribution to 
the total #-term of the theory. The next operator of dimension 6 was introduced originally 
by Weinberg |TJ]. In different classes of models violating CP-symmetry this operator may 
give an important contribution to the neutron electric dipole moment 0, [3]. 

The violation of CP-symmetry in Standard Model originates from the complexity of 
the KM matrix. To lowest, quadratic order in the weak interaction all CP-odd flavour- 
conserving amplitudes turn to zero trivially. The point is that in this approximation those 
amplitudes depend only on the moduli squared of elements of the KM matrix, so the result 
cannot contain the CP-violating phase. 

CP-odd objects may arise in the Standard Model in the fourth order in semi-week con- 
stant. However, the cancellation of EDMs of a quark and W-boson in this approximation is 
firmly established now |4], |5| . The only known non- vanishing formfactor to this approxima- 
tion is the magnetic quadrupole moment of W-boson ||. The finite EDMs can be obtained 
only after hard gluon radiative corrections are taken into account. We shall prove that in 
the absence of QCD radiative corrections the same mechanism leads to the cancellation of 
induced #-term. In contrary to the recent claim that the Weinberg operator is zero to this 
approximation |/J, we find that all operators of dim>6 acquire non- vanishing values. 

2 Schwinger operator method for calculating CP-odd 
Lagrangian 

We are going over now to the direct calculation of a few first terms of CP-odd effective 
gluonic Lagrangian in the Standard Model to three-loop approximation. The general struc- 
ture of the diagrams which could contribute to the effect in that approximation is shown 
on Fig.l, where the solid line represents a quark loop, waved lines - W-bosons. 
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The CP-odd part of the loop flavour structure reads as: 

2i5[d(c(b — s)t — t{b — s)c + t{b — s)u — u{b — s)t + u{b — s)c — c{b — s)u) 
+s(c(d — b)t — t(b — s)c + t(d — b)u — u(d — b)t + u(d — b)c — c(d — b)u) 
+b(c(s — d)t — t(s — d)c + t(s — d)u — u(s — d)t + u(s — d)c — c(s — d)u)] (2) 

For the KM matrix we use the standard parameterization of Ref. where the CP-odd 
invariant is 

5 = sin 5ciC2Cssjs 2 ss- (3) 

The letters u, d, s, c, b, t denote here the Green's functions of the corresponding quarks. 
Each product of four quark propagators allows for cyclic permutations of the kind 

udcs = dcsu = csud = sudc. 



Further considerations are based on the operator Schwinger method 0] successfully 
extended on the QCD case by Novicov, Shifman, Vainshtein and Zhakharov [JTOJ . It allows 
one to minimize the set of calculations in introducing the operator P: 




where A c (x) is the external gluonic field. Then the quark propagator taken in the back- 
ground gluonic field reads as: 

(vac\Tq a (x)q b (y)\vac) = (x, a\i{P - m) _1 |y, b) = i x ' a K P ~ m ) p2 + ig / 2 (Ga) -m 2 ^'^' 

(5) 

The external field strength appears as a result of commutation of two P: 

[P^P v ] =l G a ^=igG, v (6) 

We assume also that the field has no source and satisfies classical equations of motionF]: 

D^G^ = 0; DpGp, = 0. (7) 



The general outline of our calculation is following. Using the specific property of the 
flavour summation (|2]), namely the antisymmetrization in the masses of opposite fermionic 
lines, we rewrite the general expression for the CP-odd amplitude in KM model via some 
commutators of functions depending on P. It gives us some powers of external field and its 
derivatives in the nominator. If the explicit dimension of this part of expression is equal 
to the dimension of operator of interest, we can forget about further non-commutativity of 
Ps in other parts and put A = 0. 

The generic formula for the effective action up to some renormalization terms which 
will be discussed later looks as 



-i Tr 

flavour 



P + mi 



PP 



PP 



m 



(8) 



2 As far as we are interested in pure gluonic operators we can omit quark currents in the r.h.s. of 



3 



where the sum over quark's masses m 1; m 3 and m 2 , m 4 should be performed according 
(0). The Tr operation in this equation means the trace in colour, Lorentz and coordinate 
spaces: 

Tr[...] =Tt l+c J d 4 x(x\...\x), (9) 
T denotes the mass operator of a quark in the external field: 

f(m 2 ) = & [ l^L^ 1 + 75 *ZL 1 + 75 flu, ~ q.qjM* 

[ t) 2 J (2tt)^ 2 {q-P) 2 + ig/2{Ga)-mf v 2 q 2 — M 2 ' 1 ' 

where M is the mass of the W-boson, g w - semiweak charge. This mass operator f allows 
for the expansion in series of external field operators of increasing dimension with some 
invariant functions depending on P 2 as coefficients. It will be shown later that for the 
calculation of the Weinberg operator it is sufficient to keep in this expansion only three 
first terms of explicit dimension 0, 2 and 3.: 

fK 2 ) = (I{P,F 4 (P 2 )} + !{^(P 2 ^ (11) 

Here {...} denotes anticommutator. Some comments should be added at this point. As 



far as we do not fix the concrete view of invariant functions, two last terms in (JTT ) could 
be presented in the other form. The term of dimension 2, for instance, could be written 
as {P, {Hi(P 2 ), Ga}}. The difference between these two forms, however, is an operator 
of dimension 4 and it can be absorbed to the next term of this expansion. From the 
same reasons we do not care about the antisymmetrization in the last term in ( O) . The 
only thing which should be checked is the absence of operators {Hi(P 2 ) , {P , Ga}} and 
P> a G ^PaPpTIv This could be made from the expression (|H]) or in the framework of the 
usual perturbative expansion. 

The V - A structure of T cancels mi and 772,3 i n nominators of (§). The only prob- 
lem arises with renormalization terms which violate pure left-handedness. Now we shall 
argue that renormalization never contributes to the CP-odd effective Lagrangian to this 
approximation. 



In the absence of external field the first term in (|il|) reduces to the usual unrenormalized 
mass operator in the V-A theory: 



A=0 



P^fiP 2 ) • (12) 



The renormalization with respect to quark 1 from the left and quark 3 from the right 
introduces into the mass operator the dependence of external masses [§], |5j: 

P— 2~ /(P ) -/i3[p— 2 m i— 2 m 3^— ], (13) 

where /13 and / are expressed via the function / and masses mi, 777,3 as follows: 

ft 2n ft 2n m l h - m lh , mim 3 {fi - f 3 ) 2 2 . 

f{p)=f(P) 2 — > /i3 = 2 2 — I fi = f(P = m i)i 1 = 1,3. 14 

mf — 777,3 m i ~ m 3 
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It is clear how to adopt this scheme for the formalism of the external field. Now the first 
term of ( |TT|) should be written in the following form: 

(\{P, F(P*)} - p ^V'-mSA )!^ _£»[/>(!_ 75 ) - mi(1 - 75 ) - ros( l + 75 )], (15) 

where constants fx, / 3 and /13 are determined in (|14l). The zeroth-order term of the pertur- 
bative expansion of this expression in A reproduces fll3|); the first-order term corresponds 
to the vertex part renormalized in complying with the Ward identity. Other terms appear 
to be free of renormalization. 

Let us start our consideration of renormalization effects from the counterterms of the 
"wrong" handedness, proportional to /13. From flavour structure (Q) of the fermionic loop 
it follows, in particular, that any amplitude should be antisymmetrized in the masses mx 
and m 3 of the opposite quark's lines. The cyclic permutation allowed under the trace 
symbol simultaneously leads to the antisymmetry with respect to interchange of m 2 and 
7714. The last property automatically means the vanishing of the amplitude proportional to 
f 13(1^2} f 13(^4) ■ Thus, one of mass operators at least possesses V - A structure. Then the 
contribution of counterterms, proportional to / 13 can be easily evaluated to the form: 



Tr 



P + mx /13/A/. x x , u P + m 3 * 2 ^l + 7s 
(P(l - 75) - mi(l - 75) - m 3 (l + 75))— -V{m, 



PP 



m? 2 



PP -ml 



-A 3 TOim 3 Tr 



(p 2 + ijlGa - m\)(p 2 + i/2Ga - m§ 



-Tim 



2x 1 + 75 



• (16) 



This expression is explicitly symmetric under the permutation mi <-> m 3 and drops out 
from the answer. Therefore, only V - A structures left in both operators V and our master 
formula reduces to 



S, 



eff 



-i J2 Tr 

flavour 



P 



PP 



m\ 



Yim\) 



P 



PP -mi 



(17) 



The procedure of the antisymmetrization in mx *-> m 3 implied at this point significantly 
simplifies the set of further calculations. Indeed, the mi,m 3 -dependent part of the ampli- 
tude can be easily transformed as follows: 



P 



-r- 



p 



PP-mj PP 
P 



m% 



PP 



[T,PP] 



(mx <-> m 3 ) 
P 



P 



.PP 



m\ 



P 



PP-m\ PP -m\PP -m 



(mx <-> m 3 ) 



(PP -m\)(PP -ml) 



(m x 



m 3 l 



ml)PS 13 [t,PP]Sx3P, 



where we have introduced the operator Si 3 = (PP 



- l (PP 



mt 



~ l . It is easy to 



see that the antisymmetrization in mx <-> m 3 performed in (ff^) leads to the antisymmetry 
of the amplitude with respect to interchange of indices 2 and 4. Indeed, using the cyclic 
permutation we get: 



m 2 3 )Ti(PS 13 [T 2 ,PP}Sx 3 PT 4 



1 + 75 



-m, 



m^Tr ps 13 [r 4 ,pp]s 13 pr 2 



1 + 75 



(19) 
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The calculation of the commutator in the expression ([Tj|) is quite straightforward. For 
instance, this commutator with the rest of the first term of the mass operator can be 
calculated in the following manner: 



rl fP F(P 2 M p m iA ~ m i/3 
-{P,F{P)}-P ml _ ml 



^ = ^{A[F(P 2 ),P 2 -f^]}^ 



j({P«{F , {P lt ,D lt G aP }}} - ^{F'{P„D a D^}})^^^ + 0{dim > 5). (20) 



This formula is obtained using the relation from the general operator calculus: 

]]]}-H 

(21) 



2[F(A),B] = {F'(A), [A, B]} - ^[F"{A), [A, [A, B}]] + ±{F'"(A), [A, [A, [A, B]]]} - +... . 



The result of the commutation ( p0|) brings some important consequences: 

First, it cancels the rest of renormalization counterterms. This completes the prove of 
the statement that the renormalization never contributes to the CP-odd flavor diagonal 
amplitudes to the fourth order in semi- weak constant. 

Second, the minimal explicit dimension of the operator [f , PP] is 3 and it puts the limit 
on the number of terms in T that we have to take into account. For the calculation of the 
Weinberg operator we can restrict our considerations on the operators of explicit dimension 
3 or less and it justifies, in particular, the choice of T in the form (|TT|). 

Third, we have found the shortest way to demonstrate the absence of EDMs of the quark 
and W-boson to this approximation. There is no doubt that all considerations presented 
above can be extended on the case of the external electromagnetic field. 

Now we are in the right position for the calculation of the Weinberg operator in KM 
model. It is convenient to classify all contributions by the combination of invariant functions 
from the expansion of Ts. There are five of them: J-H, J-F, H-H, H-F, F-F. For the induced 
#-term the only possible contribution may arise from F-F. 

The simplest cases are J-H and F-H because they already possess dim=6. Therefore, 
the further non-commutativity may only influence on the effective operators of higher 
dimension. Thus, it is reasonable to make the substitutions: 

^' 13 {jp 2 — m\ )(p 2 — mj) ^ 
Then the trace over spacial variables is easily computable: 

S eff = 4g 2 Jd 4 xTr c [(D^G a p)(D u G Xf i)] x 
f d*p {m 2 ~m 2 )p a p tJi p v p x p 2 [J i (F! 2 -H 2 )-J 2 (F' i -H A )\ 



flavour 



(2tt) 4 (p 2 -mi) 2 (p 2 -ml) 



The trivial average over the direction of p produces three different field operators of the 
dimension 6. By virtue of equations of motion for the external field (0) they can be 
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transformed to the standard form of Weinberg operator. Finally, we get J-H and J-F 
contributions to the Weinberg operator in the form: 

. v f d 4 p (m 2 ~ rn^AMn - gg) ~ MFj - H A )} 
2 S LuJ W (p 2 -m 2 ) 2 ( P 2 -m 2 ) 2 1 j 



The analysis of the H-H contribution is quite transparent. The corresponding amplitude 
contains field operators of dimension 5 and 6. It means that we can neglect the non- 
commutativity of H 2 , H A with operators P and §13 because it would bring an additional 
dimension 2. Then the amplitude of interest could be transformed to the form: 

Tr({/P 2 , [0, PP]}{H 4 , 0}) = ^Tr({P 2 , [0, PP]}{H Al 6}) - {H 4 , [6, PP]}{H 2 , 6}), (25) 

where O = ig/4{P,Gcr}. It is the matter of simple exercise to check that the expression 
(f25|) vanishes identically. 

We left with the F-F and H-F groups of contributions which minimal dimension is 3. 
The straightforward calculation is quite tedious because we have to take into account Ga- 
dependence of S13, the non-commutativity of different P, etc. However, we have found a 
simple argument to show these groups do not contribute to the Weinberg operator at all. 
In the amplitude of interest (F-F case) 

Tr(PS 13 [{P, F 2 (P 2 )}, PP]S 13 P{P, P 4 (P 2 )}), (26) 

we perform the formal expansion in ig/2Ga of Fi(P 2 ) = Fi(P 2 + ig/2Gcr — ig/2Ga) around 
the "point" P 2 + igjlGo = PP. It is clear that only zeroth and first order terms of that 
expansion could be taken into account when calculating Weinberg operator. This expansion 
can be performed using another formula of the general operator calculus. Up to terms linear 
in B, function F(A + B) could be expanded in the following manner: 

2F(A + B) = 2F(A) + {F'(A), B} + ^[F"(A), [A, B]] + ... (27) 



In our case this expansion takes the form: 
F(PP - ig/2Ga) = F(PP) - ±{F'(PP),ig/2Ga} + ±[F"(PP), [PP, -ig/2Ga}} + ... (28) 

Only second term of this expansion is relevant in our consideration. Indeed, the third term 
and other denoted here as ... have the dimension higher than 3 and will not contribute 
to the Weinberg operator. First terms of the expansion of F 2 and F A drop out because 
they give a vanishing commutator with PP. So, we left with the second terms only but 
their contribution literally coincides with H-H case. Thus, the vanishing of H-F and F-F 
groups of contributions follows from the simple substitution Hi — > Hi — F[ in equation 
(p5|). Another basis of invariant functions depending on PP in the expansion of T 

tt = ({P, HP 2 )} + |{#,(P 2 ), {P, Ga}} + gJiiP^DaG^PaP^ + (29) 
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would simplify our analysis. It reduces the number of possible combinations for the calcu- 
lation of the Weinberg operator to H — H and H — J. 

We have shown that the KM-model does not induce #-term to three-loop approximation. 
To the same accuracy the Weinberg operator acquires nonvanishing contributions of the 
form flU). 



3 Weinberg operator in the KM model 

After convincing ourselves in the absence of the exact cancellation of the Weinberg operator 
in three-loop approximation, we are going to find its value. It is natural to consider all 
quark masses but m t small as compared to the W-boson one M. Together with the quark 
mass hierarchy it allows one to simplify the calculations considerably, restricting to those 
contributions to the operator of interest which are of lowest order in the light quark masses. 
Besides, it is also natural to single out the contributions with logarithms of large mass ratios, 
e.g., \og(m t /m c ), log(m&/m s ), log(M/m&) etc. 

All diagrams can be split into two types, depending on which quarks, U (u, c, t) or 
D (d, s, b), flow inside the mass operators. It is convenient to sum first of all over the 
flavours of the quarks masses of which were denoted up to now as ni\ and 777.3. For the two 
types mentioned we get respectively: 



y \ m l ~ m 3> > -m m a 

(p 2 — m 2 ) 2 (p 2 — m 2 ) 2 p 4 (p 2 — m 2 ) 2 {p 2 — m 2 ) 2 

4 2 
—m^m c 



p i {p 2 — m 2 ) 2 (p 2 — m 2 ^ 2 



(30) 



In expression ( p0|) we put m u = rrid = 0. We can determine now the characteristic momenta 
p. When quarks are arranged according to the first line of formula (0), integral fl24|) is 
infrared divergent if one neglects the masses m s and in the denominator. It means that 
the typical loop momenta contributing to the effect are p ~ m& and it cancels two powers 
of mi, m the nominator of (|30|) . In the opposite case when D-quarks are inside the mass 
operators, the typical momenta range is large: p ~ M. 

The main problem arising at this point is the calculation of invariant functions F, H 
and J. This could be done by means of usual perturbative expansion of formulae (|10"D and 
(|TT|). Using the fixed point gauge 

(z-zW(z) = (31) 

1 1 

^W) = 2 ( z ~ x )» G vn( x ) + g( 2 _ x )a( z - x) u D a G Ufl (x) + ... . 

after straightforward calculations we get the following set of equations: 

F{P > )P = P(f )P = # 1 T ^r^'r^vf - . 2 ^' (32) 



2 J (2vr) 4 (q 2 - m 2 )[(p - q) 2 - M 2 } 
H{p )G flu p fl j v — (H + F)(p )G, vP , lv = -j (27r)4(g2 _ m2)2[(p _ g)2 _ M2] ( 33 ) 
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J{ P )JJ a ^^ PuPa ^ g J (2 ^ )4 ^ _ m2 ^ p _ g y _ M 2] + 3 Vc^v&vVcnv. 

The cubic divergence of integral in the expression for F is irrelevant for us because the 
combination F'(m 2 ) — F'(m 2 ) presenting in our equations is obviously finite. 



Before taking the integrals in fl32|) - (|34| ) it is reasonable to determine the light quark 
mass dependence of F', F", H and J and sum over flavours left. Clearly, the summation 
Z)[^4(-p2 — H2) — ^2(^4 — H4)} brings an additional factor m 2 for [/-quarks, so the total 
degree of suppression for this type of diagrams is 0(m 2 m 2 c m 2 ). For the D-type of quark's 
arrangement it is essential that functions F' and F" could be expanded as follows: 



HF' 

FV, m 2 ) = F\p\ m 2 = 0) + m 2 



2^ 



+ ..., (35) 



m=0 



whereas H and J contain pieces proportional to m 2 \ogm 2 (see the Ref.0 for the details). 
In the sum over D-flavours these logarithms prevent the cancellation of terms ~ m 2 m 2 
which means that the group of diagrams with D-quarks inside mass operators contribute 
to the Weinberg operator at the same order 0(m 2 m 2 c m 2 ). Moreover, these logarithmic 
factors enhance the contribution came from D-type of quark arrangement in comparison 
with that from [/-type. Up to the last integral over p 2 the corresponding contribution to 
the Weinberg operator to double logarithmic accuracy reads as: 

^ ( 9w \ 2 mlmhn 2 , . 2/ 2 , f d 4 p (p 2 + M 2 )mf , b 2 -M 2 | 

Performing trivial integration we get the final formula for the Weinberg operator in the 
Kobayashi-Maskawa model: 

! ^^%|^logK/m 2 )log(M 2 /m 2 )/K/M 2 ), (37) 



1536tt 6 * M 4 
where the function / is 

x 2 1 / 12 10 \ x f 13 5 10 N 

I[X} ~ (x - l) 4 gX { + x - 1 + (x - l) 2 ) (x - l) 2 { + x - 1 + (x - l) 2 + (x - l) 3 

(38)' 

The Fermi constant is introduced in (|37|) according the standard notation Gf = V^g 2 / (SM 2 ] 



4 Discussion 



We have shown that the Kobayashi-Maskawa model generates CP-odd effective gluonic 
Lagrangian to three-loop approximation starting from the operator of dimension 6. The 
expression for the Weinberg operator (|37D parametrically coincides with the effective mag- 
netic quadrupole moment of H/-boson || appearing in this model in the same fourth order in 
semi- weak constant. The attempt to prove the exact cancellation of the Weinberg operator 
using the external field technic ]7) seems to be incorrect. The author of this work believes 
that antisymmetrizations of the amplitude fl8|) in mi, 1713 and 777,2, should be imposed 
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independently and both procedures increase the effective dimension. In contrary, we have 
shown they are connected and the amplitude antisymmetrized in one pair is automatically 
antisymmetric in another one. 

There is nothing surprising in the fact that the contribution to the electric dipole mo- 
ment of neutron came from ( j37D is tremendously small. For m t ~ 2M we get 

c 2 ~ 10- 27 (lGev)- 2 (39) 

Using the result of the work |Tl]] we can estimate the corresponding contribution to the 
NEDM as 10 -41 e cm. It is far beyond both current experimental limit and theoretically 



predicted NEDM came from large distances [12| as well. The extreme smallness of 



reflects not only usual parametric suppression of the effect but amazingly small numerical 
coefficient as well. 

Finally, we would like to discuss a possible value of the Weinberg operator at four- 
loop level. One additional hard gluon loop brings a factor like a s /(3ir) ~ 1CT 2 but now 
the operator of interest appears in another order in light quark masses. As a result we 
could expect the corresponding coefficient c 2 being much larger than its three-loop value 
0. To obtain an estimation for the Weinberg operator we use the approach developed 



in the work [13]. Believing that all quark masses are much smaller than the mass of W- 
boson one can use four-fermion contact limit restricting on the contributions to the effective 
Lagrangian of order M -4 . The only possible structure of diagrams is presented on the Fig. 2. 
(Dashed line here is the gluon propagator). The induced #-term appears from this graph 
in the order 0(a s G F rn 2 rn 2 ) |l3"fl . The heaviest masses enters only under logarithms here. 
The estimation for the Weinberg operator quoted in ||, 0{a s G 2 F m 2 c ), looks strange now 
because the ratio c 2 /ci is of order l/m 2 . If it is true it makes questionable the validity of the 
expansion (HD at the four-loop level. Performing the same analysis we came, however, to the 
another estimation of the effect. The typical expression for the corresponding amplitude 
before the last integral over gluon momentum k 2 looks as: 

9 9 

/Ttl Ttl 
dk 2 -^ hg(m 2 /k 2 ) \og{k 2 /m 2 c )^\og{m 2 b /k 2 ), (40) 

where the infrared divergence should be cut off at the scale of m 2 . Correspondingly, the 
four-loop contribution to the Weinberg operator is of order 0(a s G 2 F m 2 ) which is two order 
of magnitude smaller than the estimation cited above. 

I would like to thank LB. Khriplovich for helpful discussions. 
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